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Abstract 

The process of neutrinoless double electron (Oi/ECEC) capture is revisited for those 
cases where the two participating atoms are nearly degenerate in mass. The theo- 
retical framework is the formalism of an oscillation of two atoms with different total 
lepton number (and parity), one of which can be in an excited state so that mass 
degeneracy is realized. In such a case and assuming light Majorana neutrinos, the 
two atoms will be in a mixed configuration with respect to the weak interaction. A 
resonant enhancement of transitions between such pairs of atoms will occur, which 
could be detected by the subsequent electromagnetic de-excitation of the excited 
state of the daughter atom and nucleus. Available data of atomic masses, as well as 
nuclear and atomic excitations are used to select the most likely candidates for the 
resonant transitions. Assuming an effective mass for the Majorana neutrino of 1 eV, 
some half-lives are predicted to be as low as 10^^ years in the unitary limit. It is 
argued that, in order to obtain more accurate predictions for the Oz/ECEC half-lives, 
precision mass measurements of the atoms involved are necessary, which can readily 
be accomplished by today's high precision Penning traps. Further advancements 
also require a better understanding of high-lying excited states of the final nuclei 
(i.e. excitation energy, angular momentum and parity) and the calculation of the 
nuclear matrix elements. 

Keywords: neutrino mass, neutrinoless double beta decay, double electron 
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1. Introduction 

The question as to whether massive neutrinos obey a Dirac or a Majorana symme- 
try, presently constitutes one of the most important unresolved problems of particle 
physics and astrophysics. If neutrinos are Dirac particles, i.e. if neutrino and an- 
tineutrino are fundamentally different, then total lepton number L = L^. + + 
must be conserved. Contrary, if neutrinos are Majorana particles, i.e. if neutrino and 
antineutrino are identical particles, then lepton number conservation is not required 
anymore. Indeed, lepton number (LN) conservation is one of the most obscure ap- 
pearances in the Standard Model of elementary particles, since there is no known 
fundamental principle or symmetry, which would require this. 

Already in 1939 Furry [T] noticed that the exchange of neutrinos (later termed 
Majorana neutrinos) between two neutrons could lead to the production of two pro- 
tons and two electrons in the reaction 

{A,Z) ^ {A,Z + 2) + e- +e-. (1.1) 

Today, such a reaction is termed the neutrinoless double-beta (Oz//3/3) decay, and 
an observation of this reaction is still the only unambiguous way to identify the 
Majorana character of the neutrino Over the years, considerable efforts from 
experimentalists and theorists alike have been devoted to this process (for reviews 
see Ref. jsj). Although of fundamental importance, this process is unfortunately 
characterized by an excessively low rate, which poses a significant challenge to any 
experiment. 

Assuming that the light neutrino mixing mechanism provides the dominant con- 
tribution to the Oz//3/3-decay, the decay rate for a given isotope {A, Z) is simply given 
by the product of the effective Majorana neutrino mass squared |m^/3p, the known 
3-body phase-space factor, and the much less well-known nuclear matrix element 
squared, which is particular to every nuclear transition under study. The phase- 
space factor contains a dependence on the nuclear charge Z (~ Z'^), the Q- value of 
the reaction (~ Q^) and the Fermi-coupling constant (~ Gp). The main objective 
of every experimental Oz//3/3-decay search is the determination of the absolute value 
of the effective Majorana neutrino mass However, a mere observation of the 

decay would already constitute a significant advancement in neutrino physics. 

In the 3-neutrino mixing scenario, the effective Majorana neutrino mass m^^ 
takes the form 

3 

= '^U^i'n^i. (1.2) 

Here, Ugi {i = 1,2,3) are the elements of the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix, which mixes the mass eigenstates rrii in weak in- 
teraction. It contains the usual 3-neutrino mixing angles plus a CP-violating phase, 
which appears in oscillations, and two additional Majorana phases, 0i,02- 

The recent claim for an observation of the Oz//3~/3~-decay of ''^Ge with Ty^ = 
2.231 

o'.ti ^ 10^^ years implies |m^/3| ~ 0.18 — 0.30 eV by assuming the renormal- 
ized QRPA (RQRPA) nuclear matrix element and its uncertainty given in Ref. j^. 
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The goal of the upcoming GERDA experiment [6[ is to put this claim to a test by 
improving the sensitivity limit of the detection by more than an order of magnitude. 
The next generation experiments, which will be using several other candidate nuclei, 
will eventually be able to achieve this goal as well |7|. 

The value of the effective Majorana mass, as it appears in Eq. (11 ■2p . contains 
several dependencies on phases and masses. Because of experimental uncertainties, 
different mass scenarios, like the normal {m^, ^ m2 ^ mi) or inverted (m2 > nii ^ 
7713) hierarchy scenario, or degenerate {m^ ^ m2 ~ mi) or non-degenerate cases, 
can presently still be entertained, which allow a wide range of possible mass values 
for even zero in the most extreme and unfortunate situation of the normal 

hierarchy scenario jsj. Though even in that case the Oz//3/3 will still be allowed due to 
a contribution from the mass term in the neutrino propagator joj , which one usually 
neglects, its decay rate would be utterly unobservable. 

Recently, there has been an increased theoretical and experimental interest to 
another LN violating pro cess, which is the neutrinoless double electron capture 



(Oi/ECEC) [lO|, lUl, ll2|, ll3|]. In this reaction two bound electrons from the atomic 
shell are captured by two protons, thereby lowering the charge of the final nucleus 
by two units: 

{A,Z) + e,+e, ^iA,Z-2)**. (1.3) 

Here, the two asterisks denote the possibility of leaving the system in an excited 
nuclear and/or atomic state, the latter being characterized by two vacancies in the 
electron shell of the otherwise neutral atom. The energy excess given by the Q-value 
of the reaction must still be carried away by an extra photon, in order to conserve 
energy. This is unlike the 2-neutrino case, where the neutrinos can provide the 
energy balance. Thus, the reaction in Eq. (II. 3p could in principle be detected by 
monitoring the X-rays or Auger electrons emitted from excited electron shell of the 
atom, the electromagnetic decay of the excited nucleus (in case of a non-ground-state 
transition) and the extra photon, whose energy would be 

-E/ii. = Q - -£^atom - -^nucl- (1-4) 



We note that Oz/ECEC was considered by Winter [1J| already in 1955. 

The signature of the Oz/ECEC process is, therefore, different from the signature 
of the Oz//3/3 decay and would also require rather different coincident detection tech- 
niques. On the other hand, the coupling to an extra photon and/or X-ray clearly 
makes the half-life excessively long to the extent that this process has not been 
considered a valid experimental option altogether. 

The situation changes, however, if the energy difference in Eq. (II. 4p approaches 
zero and no extra photon is required. This has been discussed by Bernabeu, De 



Rujula, and Jarlskog [15[, who pointed to the possibility of a resonant enhancement 
of the Oz/ECEC decay in case of a mass degeneracy between the initial and final 
nucleus. Their best candidate case was ^^^Sn, where the Oz/ECEC double K-shell 
capture process would lead to an exited O"*" state at 1871 keV in the final nucleus 
^^^Cd. This possibility was recently excluded by a new mass difference measurement 
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performed at Jyvaskyla (AM = 1919.82(16) keV) where it was shown that 
the energy to be paid by the double K-shell vacancy would not leave enough energy 
available for the excitation of the 1871 keV state in ^^^Cd. 

The Oz/ECEC decays became a subject to a detailed theoretical treatment by 



Sujkowski and Wycech [10|, who used a perturbative approach. Their conclusion 
was that an exact energy degeneracy could make the Oz/ECEC reaction competitive 
to the Oi//3~/3~-decay. However, a case with an exact energy degeneracy could not 
be identified. 

Recently, another case for a near mass-degeneracy was found and discussed in Ref. 
ll| . Here, it was argued that the 1204 keV state in ^''Ge would be nearly degenerate 



to the ground state of the atomic nucleus ^''Se in case of a double L-capture process 
and given the experimental errors on the masses. A new mass difference measurement 
performed by Kolhinen et al. [it'] essentially confirmed the previous central mass 
difference value of AM = 1209.169(49) keV, however, with much higher precision. 
These authors therefore excluded a complete mass degeneracy with the 1204 keV 
state in ''^Ge, even in the case of a double L-capture, where the atomic energy for 



a double L-vacancy would amount to an extra 2.9 keV on top of the 1204 keV 11|. 
Prior to this, two experiments had already been performed by Barabash et al. 12l.ll3l|. 
which gave lower bounds for the half-lives of ''^Se and ^^^Sn: 

^O.ECEC(74gg) > 5.5 X 10^8 

yyECEC^ii2g^) > 9.2 X 10^9 

In this paper we present a new theoretical framework for the calculation of reso- 
nant Oi/ECEC transitions, namely the oscillation of atoms. An improved theoretical 
description of the process includes the determination of relevant matrix elements 
for the most favored cases of capture of the Si/2 and pi/2 electrons. The Oi/ECEC 
transitions without and with the spatial parity violation are considered. Further, we 
provide an updated list of the most likely resonant transitions taking new nuclear 
spectroscopic data into account and using recent accurate measurements of Q values 



for several nuclei [17|, Il6|, Il8|, |19|, |20|, |21[ . The selection of transitions is also based 
on accurate treatment of spin-tensor structures that arise in a product of the nu- 
clear matrix elements and the electron wave functions of atomic shells. The reverse 
reaction 

(A,Z)^(A,Z + 2)** + e,- + e,- (1.5) 

of a neutrinoless production of two bound electrons (Oz/EPEP) will also be discussed. 

The outline of the paper is as follows: First we discuss the mixing of atoms with 
different lepton charges. This effect leads to the oscillations of atoms. In Sect. 2, we 
discuss the relevant formalism of the oscillations. We will show that the oscillation 
of stable atoms produces a too small effect to be measured experimentally. However, 
oscillation between a stable and an excited atom can lead to a resonant enhancement 
of lepton number violating decays. 

Sect. 3 presents the estimated half-lives of the decays. In the calculations we 
use the data on the Auger and radiative widths of excited electron shells and the 
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information on the Coulomb interaction energy of two electron holes. We consider 
nuclei with arbitrary spin-parity and take into account the fact that the spin-parity 
uniquely determines a combination of upper and lower components of the relativistic 
electron wave functions entering the matrix elements associated with the capture. In 
Appendix A the procedure of averaging the electron wave functions over the nucleus 
is discussed. The transition matrix elements are derived for the JJ = 0^ and 1^ 
states of the daughter nuclei in Appendix B. The problem of calculating matrix 
elements is very complicated, and the result depends sensitively on the particular 
transition. We identify the most promising nuclei in the search for Ot^ECEC decays. 
Such nuclei will continue to be analyzed in future. In this paper, the half-lives are 
normalized to the nuclear matrix element Ai^^lO'l — )■ JJ) = 6, which is close to 
the maximum evaluated value of the matrix elements for medium-heavy nuclei. In 
Sect. 3, we also give a complete list of the most likely resonance transitions, in which 
the unitary limit of resonant enhancement gives half-lives of less than 10^^ years for 
|m/3^| = 1 eV. We argue that accurate measurements of the mass differences between 
initial and final states of the nuclei are necessary, if future experiments of Oz/ECEC 
decays with half-lives below 10^^ years were to become a possibility. Experimental 
signatures of Oz/ECEC decays are discussed in Sect. 4. 



2. Lepton number violating transitions between ground state and excited 
atoms 

If lepton number is not conserved, then the weak interaction mixing between a 
pair of neutral atoms {A, Z) and {A, Z ± 2) is a natural occurrence, which leads 
to an oscillation between these two many-body quantum systems. In the present 
description we focus on a system, in which one of the atoms (usually the daughter 
atom) is left in an excited atomic or nuclear state. In fact, for EC processes the 
daughter atomic system is always excited, as the capture process always leaves a 
vacancy in the electron shell. If the ECEC Q-value is of the order of the excitation 
of the atomic shell with two electron vacancies, one may expect a resonant-like 
transition. A few examples do exist in the nuclear chart, which have this property. 
On the other hand, if the ECEC Q-value is significantly larger than the atomic 
excitation, one may find a situation, where an excited nuclear state matches the 
available energy (i.e. Q — -Eatom = -^nuci)! allowing again a resonant-like transition 
to an excited nuclear state. The latter type of oscillations may even have a practical 
experimental signature: one or even several X-ray photons or Auger electrons from 
the de-excitation of the atomic shell being coincident with a 7-ray (or a cascade of 
7-rays) from the de-excitation of the nucleus. In fact, the detection of a coincident 7- 
ray cascade, if existent, may already be sufficient to uniquely identify the transition. 
It may be worth re-iterating that any such transition requires the neutrino to be 
of Majorana type, as there is no phase space available for the emission two extra 
neutrinos. 

The present description of a resonant enhancement of the ECEC transition will 
be done in the context of oscillations. We wish to point out, that our results are 



consistent with the results of Bernabeu et al. [I5j and Sujkowski and Wycech [10| for 
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the physically interesting case, where the frequency of oscillations is much smaller 
than the width of the excited atom. In the opposite limit, when the frequency of 
oscillations is high, the standard formulas of the time evolution of a two-level system 
are retained. 

2.1. Oscillations in arbitrary systems 

Specific features of the oscillations in the system of two atoms were discussed ear- 



lier in Ref . [22[ . Two coupled oscillators, one of which experiences friction, constitute 
the mechanical analogue of the system, which we are considering. 

Lepton number violating interactions induce transitions {A, Z) — )■ [A, Z ±2)**. 
These transitions can be described phenomenologically by 2 x 2 non-Hermitian Hamil- 
tonian matrix 

^eff = (^f: M^-fr)- (2.1) 

Here Mj and Mj are the masses of the initial and final atom. F is the decay width 
of the excited daughter atom. The off-diagonal matrix elements of i^es are complex 
conjugate. The transition potential V can always be made real by changing the phase 
of one of the states, i.e. V = V*. The diagonal matrix elements of the Hamiltonian 
are determined by strong and electromagnetic interactions, which conserve lepton 
number. The off-diagonal elements provide the mixing of the neutral atoms, and 
thereby, violate lepton number by two units as a result of the weak interaction with 
massive Majorana neutrinos. Using the Pauli matrices, the Hamiltonian can be 
written as: 

H,s = M+ + Va, + M_as, (2.2) 

where 

Mi±Mf i^ , , 

M± = ^ ^ T -f. (2.3) 

The evolution operator e''^^"^^ can be expanded over the Pauli matrices to give 



g-*iicfft = ^-^M+t l^cos(fit) - t ^ ~ sm{nt)j , (2.4) 

where n = ^V^ + Ml. 

One can see that all components of the evolution operator behave like e^''^**, 
with X± = M+ ± Q. Since the eigenfrequencies A± are complex, the norm of the 
states is not preserved in time. 

A somewhat similar form of the 2x2 Hamiltonian matrix is responsible for 



the oscillation of neutral kaons [23|. The main difference between the oscillation of 
neutral atoms and that of kaons is the mixing, which is maximum for kaons and 
exceedingly small for atoms. 

The Hamiltonian in E g. (|2.1l) also describes the effect of neutron-antineutron os- 



cillation in nuclear matter 24,|25|. In this case, Mj and Mf would be the neutron and 
antineutron masses, V the baryon number violating potential, and F the antineutron 
width related to annihilation channels (in vacuum Mj = Mf and F = 0). 

The formalism is also similar to the one used for describing oscillations and decays 



of unstable neutrinos 26 
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2.2. Oscillations of two stable atoms 

If the two atoms are stable, then F = and Q = {Mi — Mf)/2. The transition 
probabihty is determined by the off-diagonal matrix element of the evolution matrix 
(El: 

|</|e-^^-*K>P = ^sin2(fit). (2.5) 

This is just the case of oscillations of a two-level system described for instance in 
Ref. [27|. Unt ^ 1, the transition probability ~ V'^t'^ is determined by the potential 
V only. However, the exposure time of atoms in double /3-decay experiments (months 
and years) is greater than 1/Qhy many orders of magnitude. By taking the average 
over one period, we one arrives at 

In the transitions {A, Z) o {A, Z ±2), the composition of valence electron shells 
changes and, thus, the chemical properties of the substance. This circumstance can in 
principal be used for registering the oscillations of atoms. However, the potential V is 
at least 30 orders of magnitude smaller than the atomic mass difference. For a hypo- 
thetical mass difference of Mj — Mf ^ 10 keV one finds | < f\e~^^'=^^\i > p < 10~^°. 
Since degenerate ground-state masses do not exist, this scenario is purely academic, 
and we turn to systems of a stable mother and an excited daughter atom. 

2.3. Oscillations and resonant transitions between ground state and excited atoms 



According to the arguments in Ref. [22|, we assume a potential strength of ~ 
10~^^ eV, a typical decay width of F ~ 1 eV for a medium-heavy atom, and a typical 
mass difference of (Mj — Mj) ~ 1 MeV. In the lowest order in V, we obtain 

A+ = M, + AM-^Fi, (2.7) 
A_ = Mj-^F- AM+^Fi, (2.8) 

where 

Since the width Fi is small, the imaginary parts of X± in Eqs. (12. 7p and (12. 8p are 
negative, therefore, the states decay. Equation (I2.10p gives the decay rate of the 
initial atom in agreement with the Breit-Wigner formula. 

The excited atom manifests itself as a resonance in the decay amplitude. The 
most favorable conditions for the detection of the violation of lepton number conser- 
vation occur in the transitions {A, Z) — )■ {A, Z — 2)*, where the masses of the initial 
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and final states are equal (degenerate) and the decay width of the daughter atomic 
nucleus is small. 

The amplitude of finding the initial atom t seconds after its preparation in the 
same initial state is determined by the diagonal matrix element of the evolution 
operator: 

< z|e-^^-*|z >= e-^^-*^ + e-^'-\l - (2.11) 

The second term oscillates with the frequency ~ Mj and decays with the rate Fi. 
At a low frequency, the system is unable to return to the initial state and decays. 
In this case one can talk about a lepton number violating decay of the initial atom 
with a width Ti (Eq. CTID . 

The decay width Fi reaches the unitary limit 

pmax _ __ ^2.12) 

in the case of a complete degeneracy between initial and final state. From an experi- 
mental standpoint, where one would search for such lepton number violating decays, 
this would be the case of highest interest. 

3. Analysis of Oi/ECEC half-lives throughout periodic table 

The selection of atomic systems with the potentially shortest Oz/ECEC half-lives 
is based on equation f l2.10p . The equation shows that the decay rate is determined by 
three quantities: the mass difference between the initial and final states, the decay 
width of the final state, and the transition potential. 

The mass difference depends on the Q- value of the ECEC decay and the energy of 
the two electron vacancies in the final atom. In atomic physics, the electron binding 
energies can usually be calculated with an accuracy of several eV. We borrowed the 
binding energies from Ref. [28]. Noticeable corrections will, however, arise from the 
Coulomb interaction of the two holes. These calculations are carried out on the basis 
of the Dirac equation taking into account the screening effects of the nuclear charge. 

The decay width of the final atom is determined by the dipole emission rate lead- 
ing to the de-excitation of the electron shell. In the nonrelativistic approximation, 
the capture rates of two electrons from the higher shells, like L, M, or N shells, scale 
with the principal quantum numbers ni,n2 as {aZ/ni)^{aZ/n2)^, but we will see 
that in the unitarity limit of a resonant decay this strong reduction of the probabil- 
ity could possibly be compensated by smallness of the de-excitation rates. The total 
decay width of the system is given by the sum of the widths of the atomic and the 
nuclear state. In most cases the decay width of the excited nucleus is smaller than 
the one of the atomic state by at least an order of magnitude and can, therefore, be 
neglected. The process of Auger electron emission as the alternative de-excitation 



process of the atom is also taken into account following the results of Ref. [29| . The 



Auger electron emission is faster than the electromagnetic decays for low atomic Z. 
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The transition potential contains the uncertainties of the transition matrix ele- 
ments connected with the complicated structure of the nuclear excitation. In order to 
obtain numerical estimates, we factorize the Oi/ECEC matrix element on a product 
of the atomic physics factor and the nuclear matrix element. This simplification is 
justified due to weak radial dependence of the Si/2 and pi/2 electron wave functions 
inside nuclei. Further, we normalize all the Oz/ECEC nuclear matrix elements to 
the value of nuclear matrix element obtained for the ground state to ground state 
transition ^^fGd— j-^I^Sm in Ref. 3^. The contributions from the electron shell are 



determined by different combinations of the relativistic wave functions of electrons 
for the capture from different shells. We systematically examine transitions between 
all the states | — Jj| = 0, 1, 2, . . . with the parities vryvrj = ±1 for the capture of two 
electrons with orbital angular momenta Q < l\ + < 2 and the principal quantum 
numbers 1 < "^,2 < 4. 

3.1. Coulomb interaction energy of electron holes 

The binding energy of electrons in the inner atomic shells varies from 10 eV in 
light nuclei up to 100 keV in heavy nuclei. In the outer shells, the binding energy 
is a few eV, both in light and heavy nuclei. Since electrons are usually captured 
from the most inner shells, the electron binding energy gives sizeable contribution to 
the energy balance in the double electron capture. We use data of electron binding 



energies reported by Larkins [28[. Those are accurate to better than an eV for light 
nuclei and to a few tens of eV in heavy nuclei. The relevant binding energies in the 
context of ECEC are, of course, always those for the final daughter atom with Z — 2. 

In heavy elements, the electron hole interaction energies may reach values of a few 
keV (for a double K-shell vacancy), which is quite large for our problem in question. 
It is therefore essential to calculate the interaction energies of two electron holes and 
include them into the total energy equation. 

We used an approach that takes screening of the Coulomb potential by electrons 
occupying other orbitals into account. The shielding effect can be estimated from the 
known energy e of the bound electrons. In the non-relativistic theory, the effective 
charge can be found from equation 

^ = ^-^^' (3-1) 

where m is the electron mass, n is the principal quantum number. In the non- 
relativistic theory, the electron velocity v ~ aZ/n increases with the nuclear charge, 
which requires a relativistic treatment for heavy nuclei. 

The bindin g en ergies in the Coulomb field are known from the Dirac relativistic 



wave equation [3l|. Given e, the effective charge may be found from 

(3.2) 




where A = V^^^^^e^, k = -{2j + l)(j -l),j = l± 1/2, and = n - ( j + 1/2) > 0. 
Near the limit of e — )■ m, the non-relativistic formula is recovered. 
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The effective charge Z^, < Z talces into account the screening of the Coulomb po- 
tential, as well as the finite nuclear size. Given that the electron-shell wave functions 
are known, one can calculate the interaction energy of electron holes. 

We consider transitions between nuclei with good quantum numbers J'^, so the 
two-electron wave function should have good total angular momentum J, projection 
M, and parity. This can be arranged by weighting the two-electron wave function 
with the Clebsch-Gordan coefficients 

^^^f (X1,X2) = Ci'4,,„^,^/3m,(xi)vi/,^,(x2). (3.3) 

Here, \I/Q,m„(x) (a = {njl)) is the relativistic wave function of the electron in the 
Coulomb field. 

The wave function of two electrons can be written as follows: 

^ifi^u^^) = ^(^,T(xi,X2) - (-)^-''+^-^^/;^(xi,X2)). (3.4) 
The interaction energy of two electron holes can be obtained from equation 

ec= [ cixidx2^ff't(xi, X2) , , (xi, X2). (3.5) 

J |Xl — X2I 

The case of two holes with identical quantum numbers a = P and rria 7^ Tn/s requires 
special attention. The states J = 2j mod(2) are symmetric over the rria -H- 
permutations and do not exist. In the J = 2j + 1 mod(2) case the states are anti- 
symmetric over the -H- permutations, and the interaction energy should be 
divided by a factor of 2, since the superposition f l3.3p changes the overall normaliza- 
tion of the two-electron wave function, as it follows from C'-J^^ = (—l)'^~'^^C:!^^,^ ■ 

J "^IJ '"'2 ^ ' J I' •'2 J iii'V 

To simplify notations, we label the final states by indices a, 7 and the initial ones 
by The labels take values 1,2 to indicate first and second electrons. There are 
two possibilities for the final state, (a, 7) = (1, 2) and (2, 1), and two possibilities for 
the initial state, (/3,5) = (1,2) and (2, 1). 

Equation f l3.5p can be written as 



where 



and 



= - {-r^^^-' Kirsa, (3.6) 



a7/3'5 / J ^jamcj-ymj^j^mBjsms-'^ama'y'm^lSm^Sms \'^-') 

mc,m-yTnfjTns 



f e r 

/ (iXi(iX2 [^Ljxi)^^^^(xi)] \ ¥ {^2)^Sms{^2 

J |Xl — X2I L 

Further simplifications appear after the use of the expansion 

1 4:7r 

^ ^[:^nTi'^irnMY;^in2), (3.8) 



X1-X2 

im 
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Figure 1: Coulomb interaction energy of two Is electron holes versus nuclear charge Z. 



where r< = min(ri,r2), r> = max(ri,r2), and rij are unit vectors toward Xj. 
The angular integrals are calculated with the use of equation 

1 rfl]„l]Ljn)l]^„,(n)yj„(n) = ^J^Jz^^L/., (3-9) 



where Q^amSj^ are spherical spinors |27l |. 

= (-)"'*«*'"*'i/™ic;;s„ I f I I } (3.10) 

and \x\ = 2a; + 1. 

The remaining two-dimensional integral over the radial variables, 

K^f^^ = E ^^75 / ^?^^1^2^^2^-^«/.(^l)-^^^(^2), (3.11) 

can be calculated numerically or analytically using Maple. Here, 

C'X, = {-r^''^'^|W^\^ \j \ } (3.12) 

and 

K^r) = Ur)fp{r)A'^^ + g^ir)gf,,{r)A'^,p„ (3.13) 

where a = (njl), a' = {njl'), and /' = 2j — I, and similarly for /3. 

In Fig. [T] we show the Coulomb interaction energy of two K holes as a function 
of the nuclear charge. As expected, the interaction energy increases approximately 
linearly with Z. 
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More accurate estimates of the interaction energy of electron holes can be ob- 
tained on the basis of the Breit potential, which takes into account relativistic effects 
~ [y/cf. The accuracy of the present estimate, therefore, can be evaluated to be 
of order [aZjn)^. For Z = 60 the expected error in the interaction energy of two 
K electron holes is about 20%. The accuracy is, of course, better for higher-shell 
atomic excitations. 

3.2. Natural widths of excited electron shells 

Decays of excited atoms are dominated by electric dipole transitions with the 
emission of X-ray photons and/or the emission of Auger electrons. Dipole decays are 
described in the literature in detail (see, e.g., f32|). For K vacancies, characteristic 
X-ray dipole emissions are dominant in atoms with Z > 35. The dipole transition 
2p — )■ Is has a probability F = 4x 10~''Z'^ eV (see, e.g., [32|). Since for a (n+l)p — )■ ns 
decay, Ua/3 ~ 1/n'^ and dai3 ~ n^, the corresponding width decreases as Fq,^ ~ 
for transitions from higher orbits. This effect is of interest, since the unitary limit 
of the width of the lepton number violating decay is inversely proportional to the 
width of the daughter atom (see Eq. (12.121) ). 

Auger transitions of excited atoms with one electron vacancy are well studied 
theoretically and experimentally (for a review see ji^). The width of a two- hole 
state af3 is represented by 

F„/3 = F„ + F^ + F*, (3.14) 
where F* is the de-excitation width of daughter nucleus. Numerical values of the one- 



hole widths Fq, are taken from Ref. [29j. These values cover the range 10 < Z < 92 
and one-electron vacancies from K to N7 shells (1 < n < 4). Equation (I3.14p neglects 
the contributions from two-hole correlations. 

3.3. Lepton number violating potential 

We consider the Oz/ECEC process assuming the standard form of the /3-decay 
Hamiltonian 

nf'ix) = ^e{x)Y{l - ^5>e{x)Ux) + h.c, (3.15) 

where Gp = Gp cos 9c and 9c is the Cabibbo angle. The field operators of the 
electron and electron neutrino are denoted ) and Ue{x). 

The left-handed electron neutrino z/eL(x) = ^^^^^{x) is a superposition of the 
left-handed projections of Major ana neutrinos Xk with diagonal masses mk'. 

3 

UeLix) = ^UekXkLix). (3.16) 
i=l 

Here, U is PMNS neutrino mixing matrix. Majorana neutrinos are truly neutral 
particles and obey Cxk = ^kXk, where C = ^7270 is the charge conjugation matrix 
and is the phase factor. 

The strangeness conserving charged hadron current has the form 

J^(x) = pix)j''{gv - gAlMx), (3.17) 



12 



where n{x) and p{x) are the field operators of the neutron and the proton, and the 
vector and axial- vector coupling constants are gy = 1 and = 1.25. 

The potential of the Oz/ECEC capture of two electrons with the total angular 
momentum J and projection M can be written as follows 

V^,{JJ) = ^m,, (%)'-7^i^ E / rfxirfx, (3.18) 



X 



X 



+ 



g-ig-(xi-X2) 

2^^^ (2^ 

< A, Z - 2\JJxi)\n >< n\J^{x2)\A, Z > 



^.m/(Xl)C7^'^(l- 75)^/3™, (X2) J 

y 

< A,^ - 2|J^(x2)|n >< n\Jf,{xi)\A, Z > 

qo + En - Mi - Ea 



- (a ^ (3). 



Here, t/^(x) is the weak charged current in the Heisenberg representation. \A, Z > 
and \A, Z — 2 > are states of the initial and final nuclei. The sum is taken over 
all excitations of the intermediate nucleus {A, Z — 1). "^amai^) is a wave function 
of the bound electron with quantum numbers a = (riajaia), projection of the total 



angular momentum ma, and energy Ea- The factor 1/yl + takes statistics of 
the captured electrons into account: 5^/3 = 1 for the identical states and Sai3 = 
otherwise. 

In the derivation of VaplJ'^), we neglected the small neutrino masses {rrii <^ 10 
eV) in the neutrino potential, since the average exchange momentum in the process 
is large, |gj ~ 200 MeV/c. Further simplifications follows: 

i) Non-relativistic impulse approximation for the nuclear current: 

A 

J^(0,x) = Y,^-[gvg^^ + gA{cTk)ng^']Si^-^n)- (3.19) 

n=l 

ii) Closure approximation for the intermediate states: The excitation energies of 
the intermediate states En — Mj are replaced by an average value < E >^ 8 
MeV. In addition, we set Ea,i3 ~ m. The sum entering Eq. fl3.18p is then calcu- 
lated using completeness condition \n >< n\ = 1. 



iii) We restrict the calculation of the Majorana neutrino exchange potentials to the 
most favorable cases of even-even nuclei. Then the angular momentum of the 
initial nucleus is O"*" and the angular momentum of the final (possibly excited) 
nucleus J'^ must be balanced by the capture of the atomic electrons and the 
angular momentum of the atomic state. 

The potential can finally be written as 

Va^iJJ) = ^ G>/,^^v/2J7TTA1.;3(J;). (3.20) 
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Table 1: Combinations of the averaged upper and lower bi-spinor components of the electron wave 
functions entering the Oz^ECEC potential for transitions 0^ — > 0^,1^. Here, a, l3 are quantum 

numbers of electron hole states, and i^^^''(r„, r,„) and i7^^^(r„,rm) are defined in Appendix B. If 

electron holes are in the same state, Aai3 should be divided by an additional factor \/2- 



Transitions 




0+ ^ 0+ 

0+ ^ 0- 

0+ ^ 1+ 

0+ 1- 


<i^<iy(r„,rj> 
< //^y(r„,r^) > 

< (/^y (r„,r„) - (r„,r„))V4 



In the case of a capture of Si/2 and pi/2 electrons and of a favorable case for the 
nuclear transitions 0+ — )■ JJ = 0"*^, l"*^, the matrix elements M.afi{J]) are given in 
Appendix B. 

The numerical analysis of the Oz/ECEC transition is performed by factorizing the 
electron shell structure and the nuclear matrix element: 

Mc.p{J])^Ao.pM'^^{J]). (3.21) 

Here, Aap are lepton factors averaged over the nuclear volume. For low- J transitions, 
the lepton factors are given in Table 1. 

The nuclear matrix elements of 0"^ — i- 0^ transitions have the form 

M°'^(0+) = <0+||5^r-r„Mr_)[-4 + (^n-o-^)]||0+>, (3.22) 

9 A 



M'^iOj) = < 07 II J2^-r^h{r^^){K - r™) 



where 



X [^(o-„-cr„) -i(cr„ X cr„)] II 0+ >, (3.23) 

h{rnm) = -R joiqrnm) — ] ^ dq. (3.24) 

We note that the nuclear matrix element M^^{Q^) also appears in the calculation of 
the Oz//3/3-decay process [l]. 

For the transition 0"*" — )■ 1^ the lepton parts and the nuclear matrix elements 
are evaluated as discussed in Appendix A. We note that Aap vanishes for 0+ — )■ 
1^ transitions whenever the two electrons are captured from states with the same 
quantum numbers {njl) for j = 1/2. This is the reason, why e.g. the transition 
If Er(0+) Dy** (1+, E'^^^^ = 1745.72 keV) is excluded from the analysis. If 
electrons are captured from different states, e.g., two si/2 electrons from different 
shells {ua 7^ riji) or from si/2 and pi/2 states, the transition 0"*" — )■ 1^ is allowed and 
is considered in our analysis. 

The dominant combinations of upper and lower component of bi-spinors, which 
enter the lepton part of the matrix elements after the factorization, are listed in 
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Table [H The definition of functions F^\rn,rm) and H^\rn,rm) can be found in 
Appendix B. The decay rates of other transitions are estimated roughly with 

27rV2A^f, ~ ^J<f^ + 9l><f| + g}>. (3.25) 

The parity non- conservation in the weak interactions allows for instance transi- 
tions 0^ ~^ 0/ 5 1/ accompanied by capture of two S'-wave or two P-wave electrons. 

3.4- Likely resonant OuECEC transitions 

We have considered all the nuclei and their excited states registered in the 



database of the Brookhaven National Laboratory [33|] in August 2010, as well as 
all the combinatorial possibilities associated with the capture of two electrons. The 
selection criteria are as follows: 

i) The excitation energies are usually known with precision much higher than the 
atomic ground-state masses. We selected those pairs, where degeneracy occurs 
within the bounds given by a three standard deviation error of the ground-state 
mass measurements. 

ii) The unitary limit for the normalized half-life is less than 10^^ years. 
The half-lives are calculated using the formula 

ln2 

Ti/2 = -f^, (3.26) 
J- 1 

where Fi is given by Eq. (12.101) . 

Tables [2], [3], and H] present the results of the selection over stable parent isotopes. 
We show the natural abundances NA, the spin-parity of the final nucleus JJ, the 
excitation energy of the final nucleus E^^^i = MXz-2 ~ Ma_,z-2, and the total mass 
difference M^*^_2 — Ma,z- The two errors indicate the errors of the ground-state 
mass measurements. Shown are as well the quantum numbers of the two hole states 
a and /3 in the electron shell, the energy of the holes (not including the electron 
rest mass), the Coulomb interaction energies of the holes ec, and the decay widths 
Fq,^. The last two columns show the minimum and maximum normalized Oz/ECEC 
half-lives. 

Tables [21 [3l and H] list all, but no more than 5 transitions with the lowest quantum 
numbers of electron holes for each pair of the elements. If the spin is not fully 
determined, we took its lowest suggested value. 

Rigorous calculations of the nuclear matrix elements (NME) based on the struc- 
ture of nuclear states have not yet been performed. The objective here is to first select 
promising pairs of nuclei on the basis of rough estimates of the matrix elements, as 
these won't significantly change the global picture. The half-lives are normalized to 
the nuclear matrix element of Ai^^{0'j) = 6, which roughly corresponds to the max- 
imum evaluated value of NMEs for medium- heavy nuclei |3] . Transitions to excited 



states are suppressed due to dissimilarity of the nuclear wave functions |3J, [35 
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Tabic 2: Likely resonant Oi^ECEC transitions. The list contains only those initial and final atoms, which are potentially degenerate at the level of 
a 3(7 experimental error, and for which the unitary limit of the half-lives is below 10^^ years. The first column shows the natural abundances 

(NA in %) of the parent nuclei. In column 3 the spin and parity of the excited final nuclei (or their suggested values in parentheses or, if unknown, 
their assumed values in square brackets) are listed. Column 4 shows the excitation energies of the final nuclei with their present experimental errors. 
Column 5 lists the total mass differences, including the hole energies of the final nuclei, where the first error indicates the experimental uncertainty of 
the parent ground-state mass and the second the one of the daughter. The quantum numbers of the electron holes in the next two columns are {n2jl), 
which arc the principal quantum number, twice the total angular momentum, and the orbital angular momentum. Columns 8, 9, and 10 list the hole 
energies and the Coulomb interaction energies. Column 11 shows the widths of the excited electron shells. The last two columns show the minimum 
and maximum half- lives of the Oi/ECEC transitions (in years). Masses, energies and widths are given in keV. 



NA 


Transition 




^A,Z-2 - Ma,Z-2 




(n2jl)a 


{n2jl)p 










'T^min 
-'1/2 


T^max 
^1/2 


5.52% 


9|Ru-)-|6Mo** 


0+ 


2742 ± 1 


24.1± 7.9± 1.9 
23.7± 7.9± 1.9 


310 
410 


410 
410 


0.50 
0.06 


0.06 
0.06 


0.02 
0.01 


9.5 X 10"^ 
6.4 X 10--^ 


3 X 1026 

1 X 1027 


9 X 1033 
6 X 1034 


1.25% 




[0+] 


2737 ± 1 


16.5± 5.9± 4.1 
-4.8± 5.9± 4.1 
-5.1± 5.9± 4.1 
-7.9± 5.9± 4.1 
-8.5± 5.9± 4.1 


110 
110 
110 
110 
110 


110 
210 
211 
310 
410 


24.35 
24.35 
24.35 
24.35 
24.35 


24.35 
3.60 
3.33 
0.67 
0.09 


0.74 
0.23 
0.21 
0.07 
0.02 


1.3 X 10"^ 

1.0 X 10"^ 
8.5 X 10-3 

1.4 X 10-2 

1.1 X 10-2 


3 X 1023 
9 X 1023 
1 X 1026 

4 X 1024 

7 X 1024 


2 X 103" 

3 X 1030 

5 X 1032 

9 X 1030 

3 X 103i 


0.095% 


1 2 4 v £1 _i. 1 2 4 'T'„ * * 
54 


[0+] 


2853.2 ±0.6 


-1.2± 1.8± 1.5 
-1.6± 1.8± 1.5 
-5.2± 1.8± 1.5 
-5.4± 1.8± 1.5 
-6.1± 1.8± 1.5 


210 
210 
210 
210 

210 


210 
211 
310 
311 

410 


4.94 
4.94 
4.94 
4.94 

4.94 


4.94 
4.61 
1.01 
0.87 

0.17 


0.16 
0.16 
0.08 
0.06 
0.02 


4.4 X 10-3 
5.0 X 10^3 
1.2 X 10-2 
5.4 X 10-3 

4.6 X 10^3 


2 X 1024 

2 X 1026 

9 X 1024 
5 X 1026 

8 X 102-* 


3 X 1030 

2 X 1032 

8 X 1030 

2 X 1033 

7 X 103i 


0.185% 




0+ 


2315.32 ±0.07 


-27.5±13.3± 0.4 


110 


110 


37.44 


37.44 


0.93 


2.6 X 10^2 


1 X 1023 


7 X 1023 


0.185% 




[0+] 


2349.5 ±0.5 


6.7±13.3± 0.4 
-25.4±13.3± 0.4 
-25.8±13.3± 0.4 
-30.3±13.3± 0.4 
-30.5±13.3± 0.4 


110 
110 
110 
110 
110 


110 
210 
211 
310 
311 


37.44 
37.44 
37.44 
37.44 
37.44 


37.44 
5.99 
5.62 
1.29 
1.14 


0.93 
0.30 
0.28 
0.09 
0.08 


2.6 X 10^2 

1.5 X 10-2 

1.6 X 10-2 
2.4 X 10-2 

1.7 X 10-2 


1 X 1023 
3 X 1023 
3 X 1025 
1 X 1024 
1 X 102'5 


2 X 1023 

4 X 103" 

4 X 1032 

1 X 1031 

2 X 1033 


0.185% 




(l+,2+) 


2392.1 ±0.6 


17.1±13.3± 0.4 
16.8±13.3± 0.4 
12.3±13.3± 0.4 
12.1±13.3± 0.4 
11.2±13.3± 0.4 


110 
110 
110 
110 
110 


210 
211 
310 
311 
410 


37.44 
37.44 
37.44 
37.44 
37.44 


5.99 
5.62 
1.29 
1.14 
0.25 


0.21 
0.29 
0.07 
0.08 
0.03 


1.5 X 10-2 

1.6 X 10-2 
2.4 X 10-2 

1.7 X 10-2 
1.6 X 10-2 


8 X 1022 

2 X 102-* 

4 X 1023 
7 X 1024 
7 X 1023 


6 X 1023 

2 X 103i 
9 X 1029 

3 X 103i 
3 X 1030 



Table 3: Continued from Table El 



NA 


Transition 


JTV 




MT,z-2 - Ma,z 


{n2jl)a 












^1/2 


-'1/2 


0.185% 




(l+,2+) 


2399.87 ±0.05 


24.9±13.3± 0.4 
24.6±13.3± 0.4 
20.0±13.3± 0.4 
19.9±13.3± 0.4 
18.9±13.3± 0.4 


110 
110 
110 
110 
110 


210 
211 
310 
311 
410 


37.44 
37.44 
37.44 
37.44 
37.44 


5.99 
5.62 
1.29 
1.14 
0.25 


0.21 
0.29 
0.07 
0.08 
0.03 


1.5 X 10"^ 

1.6 X 10-2 
2.4 X 10-2 

1.7 X 10-2 
1.6 X 10-2 


8 X 1022 
2 X 1024 
4 X 1023 
7 X 1024 
7 X 1023 


1 X 10^" 
3 X 10^1 

2 X 10^" 
5 X 10^1 
5 X 10^" 


0.20% 




0+ 





-0.8± 2.5± 2.5 
-1.3± 2.5± 2.5 
-7.1± 2.5± 2.5 
-7.3± 2.5± 2.5 
-8.5± 2.5± 2.5 


110 
110 
110 
110 
110 


210 
211 
310 
311 
410 


46.83 
46.83 
46.83 
46.83 
46.83 


7.74 
7.31 
1.72 
1.54 
0.35 


0.34 
0.32 
0.11 
0.10 
0.04 


2.3 X 10-2 

2.3 X 10-2 
3.2 X 10-2 
2.5 X 10-2 

2.4 X 10-2 


2 X 1023 

2 X 1025 

7 X 1023 

5 X 1025 

1 X 1024 


8 X 1026 
2 X 1029 

1 X 1029 

2 X 10^1 
6 X 1029 


0.06% 






1946.375±0.006 


-7.5± 6.6± 2.5 
-14.0± 6.6± 2.5 
-15. 5± 6.6± 2.5 


110 
110 
110 


211 
311 
411 


50.24 
50.24 
50.24 


7.93 
1.69 
0.29 


0.34 
0.10 
0.03 


2.6 X 10-2 
2.8 X 10-2 
2.8 X 10-2 


8 X 1023 

2 X 1024 

5 X 1024 


5 X 1029 
3 X 10^0 
8 X 10^0 


0.06% 




0- 


1952.385±0.007 


-1.1± 6.6± 2.5 
-1.5± 6.6± 2.5 
-7.8± 6.6± 2.5 
-8.0± 6.6± 2.5 
-9.4± 6.6± 2.5 


110 
110 
110 
110 
110 


210 
211 
310 
311 
410 


50.24 
50.24 
50.24 
50.24 
50.24 


8.38 
7.93 
1.88 
1.69 
0.38 


0.35 
0.34 
0.11 
0.10 
0.04 


2.6 X 10-2 

2.6 X 10-2 
3.5 X 10-2 
2.8 X 10-2 

2.7 X 10-2 


2 X 1024 

4 X 1025 
9 X 1024 

1 X 1026 

2 X 1025 


6 X 1029 
1 X 10^1 
3 X 10^" 

7 X 10^1 
1 X 10^1 


0.06% 


efDy^ef Gd" 




1962.037±0.012 


8.1± 6.6± 2.5 
1.6± 6.6± 2.5 
0.2± 6.6± 2.5 


110 
110 
110 


211 
311 
411 


50.24 
50.24 
50.24 


7.93 
1.69 
0.29 


0.34 
0.10 
0.03 


2.6 X 10-2 
2.8 X 10-2 
2.8 X 10-2 


8 X 1023 
2 X 1024 

5 X 1024 


5 X 1029 

6 X 1029 
1 X 10^0 


0.06% 


ei'Dy^ef Gd" 


1+ 


1965.950±0.004 


12.4± 6.6± 2.5 
12.0± 6.6± 2.5 
5.7± 6.6± 2.5 
5.6± 6.6± 2.5 
4.2± 6.6± 2.5 


110 
110 
110 
110 
110 


210 
211 
310 
311 
410 


50.24 
50.24 
50.24 
50.24 
50.24 


8.38 
7.93 
1.88 
1.69 
0.38 


0.26 
0.34 
0.09 
0.10 
0.03 


2.6 X 10-2 

2.6 X 10-2 
3.5 X 10-2 
2.8 X 10-2 

2.7 X 10-2 


4 X 1022 
8 X 1023 

2 X 1023 

2 X 1024 

3 X 1023 


4 X 102** 

8 X 1029 

4 X 1028 

9 X 1029 
1 X 1029 


0.06% 


efDy^ef Gd" 


[0+] 


1970.2 ±0.8 


16.7± 6.6± 2.5 
16.3± 6.6± 2.5 
10.0± 6.6± 2.5 
9.8± 6.6± 2.5 
8.4± 6.6± 2.5 


110 
110 
110 
110 

no 


210 
211 
310 
311 
410 


50.24 
50.24 
50.24 
50.24 
50.24 


8.38 
7.93 
1.88 
1.69 
0.38 


0.35 
0.34 
0.11 
0.10 
0.04 


2.6 X 10-2 

2.6 X 10-2 
3.5 X 10-2 
2.8 X 10-2 

2.7 X 10-2 


1 X 1023 
1 X 1025 
6 X 1023 
4 X 1025 
1 X 1024 


2 X 1029 

2 X 10^1 

3 X 1029 

3 X 10^1 

7 X 1029 



Table 4: Continued from Tabled 



NA 


Transition 




MXz_2-Ma,z-2 




{n2jl)a 


{n2jl)fi 










^1/2 


-'1/2 


0.06% 




0+ 


1988.5 ±0.2 


-7.0± 6.6± 2.5 
-7.4± 6.6± 2.5 


210 
210 


210 
211 


8.38 
8.38 


8.38 
7.93 


0.21 
0.22 


7.6 X IQ-^ 

7.7 X 10-3 


5 X 1023 
2 X 1025 


3 X 103f 

2 X 1032 










-13.6± 6.6± 2.5 
-13. 8± 6.6± 2.5 
-15. 1± 6.6± 2.5 


210 
210 
210 


310 
311 
410 


8.38 
8.38 
8.38 


1.88 
1.69 
0.38 


0.10 
0.09 
0.03 


1.7 X 10-2 
9.4 X 10-"^ 
8.7 X 10-3 


2 X 102-1 
8 X 102''' 

2 X 1024 


6 X 1030 

8 X 1032 

3 X 103i 


0.06% 




1+ 


2026.664±0.006 


17.9± 6.6± 2.5 
16. 5± 6.6± 2.5 
16.4± 6.6± 2.5 
16.3± 6.6± 2.5 
16.2± 6.6± 2.5 


310 
310 
310 
311 
311 


311 
410 
411 
410 
411 


1.88 
1.88 
1.88 
1.69 
1.69 


1.69 
0.38 
0.29 
0.38 
0.29 


0.07 
0.02 
0.03 
0.03 
0.02 


1.8 X 10-2 
1.8 X 10-2 
1.8 X 10-2 

1.0 X 10-2 

1.1 X 10^2 


3 X 1025 

4 X 102'1 

7 X 1025 
4 X 1025 

7 X 102*5 


1 X 1032 

2 X 103i 

3 X 1032 
5 X 1032 
8 X 1033 


1.56% 




0+ 





-5.5± 3.1± 2.5 
-5.9± 3.1± 2.5 
-12.6± 3.1± 2.5 
-12.8± 3.1± 2.5 
-6.4± 3.1± 2.5 


210 
210 
210 
210 
211 


210 
211 
310 
311 
211 


9.05 
9.05 
9.05 
9.05 
8.58 


9.05 
8.58 
2.05 
1.84 
8.58 


0.22 
0.23 
0.11 
0.09 
0.27 


8.6 X 10^3 
8.3 X 10-3 
1.8 X 10-2 
1.0 X 10-2 
8.0 X 10^3 


4 X 1023 

2 X 1025 

1 X 102-1 
6 X 1025 
1 X 102fi 


1 X 103f 

5 X 103i 

3 X 1030 

4 X 1032 

5 X 1032 


0.14% 


i«8Yb^i|«Er" 


1- 


1358.898±0.005 


3.9± 4.4± 2.5 
-3.6± 4.4± 2.5 
-5.3± 4.4± 2.5 


110 
110 
110 


211 
311 
411 


57.49 
57.49 
57.49 


9.26 
2.01 
0.37 


0.38 
0.11 
0.04 


3.3 X 10-2 
3.5 X 10-2 

3.4 X 10-2 


4 X 1023 
1 X 1024 

3 X 102-1 


7 X 1028 

2 X 1029 
5 X 1029 


0.14% 


leSYb^lfEr" 


0+ 


1422.10 ±0.03 


12.0± 4.4± 2.5 
11. 8± 4.4± 2.5 
11. 5± 4.4± 2.5 
11.4± 4.4± 2.5 
9.7± 4.4± 2.5 


210 
210 
211 
211 
211 


310 
311 
310 
311 
410 


9.75 
9.75 
9.26 
9.26 
9.26 


2.21 
2.01 
2.21 
2.01 
0.45 


0.11 
0.10 
0.10 
0.13 
0.04 


1.8 X 10-2 
1.1 X 10-2 
1.8 X 10-2 
1.1 X 10-2 
1.0 X 10-2 


1 X 102-1 
4 X 1025 

8 X 102^ 

2 X 1026 

1 X 1026 


2 X 1030 

2 X 1032 
2 X 1032 
8 X 1032 
5 X 1032 


0.13% 


If W^ifOHf* 


0+ 





-12.0± 3.9± 2.1 


110 


110 


65.35 


65.35 


1.26 


7.2 X 10-2 


3 X 1022 


4 X 102'^ 


0.02% 




(0) + 


1322.152±0.022 


11. 3± 1.3± 0.9 


110 


110 


69.53 


69.53 


1.31 


8.0 X 10-2 


7 X 1026 


2 X 102'^ 


0.014% 


i|Opt^i|OOs" 


(1,2,3) + 


1382.4±0.2 


15.6± 5.7± 1.5 
15.4± 5.7± 1.5 
13.2± 5.7± 1.5 
3.2± 5.7± 1.5 
0.8± 5.7± 1.5 


210 
210 
210 
310 
410 


310 
311 
410 
410 
410 


12.97 
12.97 
12.97 
3.05 
0.65 


3.05 
2.79 
0.65 
0.65 
0.65 


0.13 
0.12 
0.05 
0.04 
0.03 


2.2 X 10-2 
1.6 X 10-2 
1.5 X 10-2 
2.2 X 10-2 
1.5 X 10-2 


3 X 1023 
1 X 1025 
5 X 1023 

3 X 102-1 
9 X 102-1 


8 X 1023 
6 X 103i 
2 X 1030 

1 X 103f 

5 X 1030 



The half-lives are further normalized to an effective neutrino mass of |?7i^/3| = 1 eV. 
The normalized half-lives listed in Tables 121 121 and H] are then given as 



2^1/2 = ^1/2 



leV 



(3.27) 



The unitary limit, i.e. the minimum value of T1/2 (denoted as T™™) is given for a 
full mass degeneracy of initial and final atoms. The maximum value is obtained 
by substituting in Eq. fl2.10p the mass difference squared 

AM^ = {MXz-2 - MA,zf + AM,2,p„ (3.28) 

where the last term accounts for the experimental errors of the masses and the nuclear 
excited state: 

AM,%, = 5Ml,_, + 5Ml, + 5Ell,. (3.29) 

The atomic mass differences are occasionally known with higher precision than the 
individual masses. The normalized half-lives, however, should stay between the two 
bounds given by f^'^ and f^'^f . 

There are some other Oz/ECEC transitions, which may deserve an extra comment: 

• 3lKr(0+) ^ISe** 

This transition has a decay ECEC Q- value of 2867.5 keV and appears in the list 
of likely resonant transitions in Ref . . The final nucleus has an excited state 
at 2864 keV, however, with unknown spin and parity. The state, however, 7 
decays by 100% into the 3"*" state at 1054 keV state, which essentially excludes 
a spin J = 0, 1 assignment for the 2864 keV state. This transition is excluded 
from our list. 

• ^^Ru(0+) ^^6mo** 

A mass degeneracy of initial and final atoms could occur for this transition 
by considering the excitation of final nucleus with 2712.68 keV. The angular 
momentum and parity of this excited state are not known yet. However, it 
decays to state with = 5^/6^, which excludes a spin J = 0, 1 assignment 
for this state. Therefore, this transition is excluded from the list as well. 



i06Cd(0+) ^ll^'Pd** 

This transition might be a resonant transition in the case the palladium iso- 
tope remains in an excited nuclear state at 2717.59 keV. However, the state 7 
decays by 100% into the 3"*" state at 1557.68 keV state, which again excludes 
a possibility of J = 0, 1 for this state. 

^i2Sn(0+) ^ll'Cd**{0+) 

The transition through the excited level at 1871 keV (O"*") has been widely 
discussed in the literature [IB] . As already noted earlier, the masses and thereby 
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the Q- value for the ground-state ECEC decay have been re-measured with high 



accuracy to Q = 1919.82 keV [16|. A double K-shell excitation of the final atom 
requires about 55 keV, thereby leaving not enough energy to excite the 1871 keV 
state. This transition is therefore not considered anymore in Tables [21 [3l and|H 

The excitation level of 2544.43 ±0.08 keV with unknown spin-parity gives a 
quite low half-life assuming = 0^ and a double capture of K electrons. If 
higher electron levels are involved, the degeneracy gets broken. The excited 
IfXe* level decays by 100 units into the 2+ state at 536 keV and by 10 units 
to the ground state. The strong ground state mode excludes a 0^ assignment. 
Possible assignments are J > 1. Two K electrons have J = 0, so such options 
are forbidden by the conservation of angular momentum. 

There is a possibility of a resonant enhancement of the Oz/ECEC of ^^°Ba by 
considering a nuclear excitation of ^^°Xe at 2608.43 keV or 2622.32 keV. The 
angular momentum and parity of these states are not clarified yet. Unfortu- 
nately, a favored possibility of 0"*" excited state is excluded, because the 7-decay 
of these levels feed into a 4+ state. These transitions have not been considered 
in Tables m-H 

• ^rEr(0+) Dy-(1+) 

The decay of Er (ECEC Q- value 1845 keV) with the excitation of the 1745.72 
(JJ = 1"*") keV state in ge^Dy is a good candidate as far as energy matching 
is concerned. A double K-shell capture requires about 110 keV, however, the 
transition to a I''" final state is strictly Pauli-forbidden. 

• If Os-^ifW**(0)+ 

In the transition to the 1322 keV excited state of 74^W, the mass difference 
exceeds three standard deviations. Accordingly, we do not expect a complete 
degeneracy. However, the minimum half-life estimated for three standard de- 
viations in the direction of a smaller mass difference is less than 10^^ years. 
This decay is included in the tables. Shown in place of the unitary limit is the 
minimum half-life found in this way. 

All systems in Tables |2l [3l and H] start out from stable parent nuclei. We also 
investigated unstable and radioactive candidate nuclei, in particular those with rather 
long half-lives. However, none of them was found to be of practical use. Even if 
those nuclei could be produced in reasonably large quantities, activity levels would 
be prohibitive. 

3. 5. Likely resonant OuEPEP transitions 

The selection criteria have also been applied to the neutrinoless double electron 
production (EPEP), where the two electrons are placed into a bound state above the 
occupied electron shell of the final atom (see Eq. (11. 5p ). and where the simultaneous 
nuclear excitation provides the mass degeneracy. Clearly, this process is expected to 
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be rather unlikely, as it requires that a Q- value be extremely fine-tuned to a nuclear 
excitation. If such a case existed, the two electrons would be placed into any of the 
upper most non-occupied electron shells of the final atom, whereby the width of the 
resonant decay would then be controlled by the width of the excited nucleus. The 
atomic de-excitation width would be comparatively small, or even zero in the case 
the electrons occupied the atomic ground state. 

As an instructive test, we analyze the Oz/EPEP transition of gQ^Nd(0"^) (isotopic 
abundance: 5.8%) to the 1920.97 keV excited 0+ state in g2^Sm. In this case the 
difference of the atomic masses is 2,+2 ~ ^a,z = 7.2 ± 2.8 ± 2.4 keV. The ground- 
state configuration of the 62^Sm daughter atom is 4/^6s^, whereas the one of go^Nd 
is 4/^6s^. The Oi^EPEP transition to the 4/^6s^ ground state of 62'^Sm is disallowed 
by angular momentum, as it requires the production of two electrons with / = 3. The 
most favored Oi/EPEP transition would be an atomic excitation of 62*Sm, whereby 
the two electrons were placed into the 7s shell. The atomic width is small, and 
we estimate that the unitary limit assuming the nuclear half-life of 0.1 ps. The 
calculated shortest half-life would be above 10^'' y. 

3. 6. Data analysis of candidate transitions 

Tables [21 [3l and H] also give a list of decays in which the unitary limit for the half- 
life of the atom turns out to be low, because the low probability of capture of electrons 
from higher orbits is compensated by a low probability of de-excitation of the electron 
shell. The probability of finding degeneracy with an accuracy comparable to the de- 
excitation width, of course, decreases with the de-excitation width. The number 
of such transitions, however, is sufficiently large, so it makes sense to evaluate the 
chances of detecting degeneracy with improving the accuracy of the measurements, 
taking into account all the states. We restricted ourselves to the capture of electrons 
from the orbits with principal quantum number no higher than 4. Higher states may 
also participate in the decays, however, we have not considered them, because data 
on the de-excitation widths of electron shells with principal quantum number higher 
than 4 are not available. 

The objective of the analysis of this subsection is to estimate a priori probability 
of finding in future experiments an atom with a half-life of less than some fixed value. 

The analysis is made by considering all transitions listed in Tables |2l [3l and HJ 
as well as the transitions associated with the capture of electrons from higher orbits, 
that passed our filters, but were not included in Tables [H El and IH as long as we list 
not more than 5 transitions for each pair of the nuclei, associated with the electron 
capture from lowest orbits. 

In the spirit of the Bayesian method, the mass difference between initial and 
final atoms is regarded as a random variable. We assume further that this value, 
?7, is normally distributed around the current experimental value M^*^_2 — M^^z 
with a dispersion determined by the current experimental error AMexpt- Future 
measurement of the mass difference between pairs of atoms with an accuracy of A/x 
is consistent with the hypothesis of the degeneracy of atoms with a probability 



w(A/i) = 




Am v27rAMe 



expt 



exp(— 



2AM2 



(3.30) 
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Suppose a measurement gives 77 = with the accuracy of A^u ~ Ta/s- For the 
normal distribution of 7], in the hkely resonant case {M'^^_2 — M^^zY ~ ^^expf 
and for a small Tap-, the probability of finding the complete degeneracy and thereby 
of getting the unitary limit equals w(A/i) ~ Fa^/AMexpt, in agreement with naive 
expectations. 

The half-live T1/2 of a particular transition may be written as 

Ti/2 = Ti7- , (3.31) 

where AM is defined in Eq. f l3.28p . To have a decay time less than T1/2, it is sufficient 
to claim experimentally 77 = with an accuracy better than 

A,<(f,,.) = I|£ (^ll - 1 j . (3.32) 

The value of Ayu(Ti/2) is the maximum error, at which the decay time is shorter than 

^1/2- 

The number of transitions, n^, to be found with half- lives below T1/2 can be 
estimated by summing up the Gaussian probabilities 

nT = ^w(A/i(Ti/2)) 



The sum runs over the transitions with T™™ < T'i/2; listed in Tables |2l [3l and|H 

The results are shown in Fig. [2l There are 2 chances out of 100 to find the 
transition with a normalized half-life below 10^^ years, 13 chances out of 100 to find 
the transition with a normalized half-life below 10^^ years, 45 chances out of 100 to 
find the transition with a normalized half-life below 10^® years. We also expect ~ 2 
transitions with T1/2 < 10^'' years. 

Under the specified conditions, the chance of finding after precise mass measure- 
ments a Oz/ECEC decay more sensitive to the Majorana neutrino mass than the 
benchmark ^^Ge neutrinoless (3~(3~ decay is close to ~ 1 : 25. The experimental 
search for the Oz/ECEC decay is thus expected to be more complicated. The back- 
ground conditions in the Oz/ECEC decay are, however, very favorable (see Sect. 4). 

Here, all the transitions are treated as statistically independent. For most of the 
transitions in Tables [21 [3], and HI this condition is certainly not satisfied, because 
for the same atom, together with the capture of electrons from the lower shells, one 
generally has captures of electrons from the high-lying levels. One precise measure- 
ment of the atomic mass difference either excludes the entire group, or allows all (or 
nearly all) transitions. Breit-Wigner resonances corresponding to the different elec- 
tron states have a typical width of less than ~ 10 eV. In our case, such resonances do 



V- / 1 MXz-2-Ma,z + A^(Ti/2) 

^ V 2^^ ^ y2AMe: 



^cxpt 



1 M*/y 2 - z - Au(Ti/2) , 
--erf( ^ ' . (3.33) 
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Figure 2: Expected number of transitions with a normalized half- life below T1/2, provided ac- 
curate measurements are made to determine the mass difference between pairs of atoms listed in 
Tables El El andlH 



not overlap and are summed up additively. An increase in the decay probability is 
equivalent to registering an excess of events. The number of events calculated from 
Eq. (13.331) takes into account this effect. Although the list of transitions splits into 
the correlated groups, the statistically independent count of the expected decays is 
justified. While the correlations do not influence tit, they obviously increase the 
dispersion of n^. 

The accuracy of 10 eV in the measurement of atomic masses will be achievable 



in the near future [18j. The electron binding energy depends on the local physical 
and chemical environment. An interesting question is whether and how to manage 
the atomic structure in such a way as to artificially implement the degeneracy of the 
atoms and create conditions for the resonant enhancement. 

The effects of finite density and temperature on the profile of the atomic spectral 
lines were studied theoretically and experimentally starting from the 1930's. To give 
an idea of the magnitude of the effect, we give some numbers. In ionized atoms of He 
I, photon energy emitted in the transition between the levels 3p and 2s at a density 
of 10^®cm~^ and temperature T = 5000 K changes by about 10 meV (see, e.g., [3^ . 
Sect. 39). In metals, when the temperature varies from 100 K to room temperature, 
levels Is in Li, 4p in Rb get shifted by several tens of meV [40;]. These values 
are 2-3 orders of magnitude smaller than the width of electron holes in the atoms, 
so the influence of the environment in such a case is not sufficient and, apparently, 
is not of interest. 

The chemical composition of matter also may influence the binding energies. In 
comparison with the free atoms, the energy of atomic electrons in metals gets shifted 



by a few eV [41|. The total effect on the mass difference of atoms connected by 



double-electron capture can reach ~ Z eV, which exceeds the typical natural width 
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~ 10 eV of the electron holes in medium-heavy and heavy atoms. The effect of 
chemical composition deserves more detailed analysis. 



4. Experimental signatures of Oi/ECEC decay 

Double beta decay of any variant is a process with a notoriously small 

probability, where typical decay times are at least of the order 10^^ y in the most 
fortuitous cases. These values pertain to the ordinary cases, where 2v are involved. 
For the more relevant neutrinoless decays, half-lives are expected to be even longer 
by at least another 3 to 4 orders of magnitude owing to the fact that the rate scales 
with the square of the mass of the Majorana neutrino. To detect these extremely rare 
processes in a real experiment, requires an enormous effort in background reduction. 
There are two typical and rather different sources of backgrounds one has to deal 
with, cosmic ray interactions (including cosmogenically produced radio-isotopes) and 
the ambient radioactivity. Even in underground laboratories like the ones at Gran 
Sasso or at Modane, which provide a natural 3100 meter, resp. 4800 meter water 
equivalent shielding against cosmic rays, an experimental setup for decay mea- 
surements still requires a massive active and passive shielding before any attempt 
can be made to venture into life-times of the order of 10^^ y and above p, |36|, [37 . 
Typical probe masses are then still several tens or hundreds of kilograms. 

Double beta decay experiments so far have focused almost entirely on the 
decay variant and little attention has been devoted to a possible decay from the 

direction. Indeed, the /3~/3~-decay comes with a rather simple and compara- 
tively easy to detect signature: in the 2i//3~/3~-decay the summed energy spectrum 
of the two electrons carries the signature of a 4-body decay, whereas a single, mono- 
energetic peak at the endpoint energy signals the neutrinoless decay. However, be- 
cause the summed electron energy of the 2i//3~/3~-decay extends all the way to the 
endpoint region, any background in this region becomes a serious issue. The spectral 
resolution of the experiment is then most critical parameter, which determines the 
sensitivity of the experiment. 

The signatures of a decay from the /3~^ direction are rather different. We will 
not enter into a detailed discussion about (3~^(3~^ or /3+EC decay, as there is an 
extra energy penalty of 2me(? to be paid for each Z?"*" production on top of a rather 
disadvantageous (i^ phase-space factor and final-state Coulomb repulsion. 

We focus on the ECEC process instead. In comparison with the /3~/3~-decay, 
prospects for measuring the ECEC-decay looked at the beginning rather pessimistic. 
Generally low isotopic abundances add to the difficulties. At present, there is only 
one experiment, which has advanced quite significantly. This is the TGV experi- 



ment [38|] in the Modane underground laboratory looking for the ECEC decay of 
^"^^Cd through the identification of ~ 20 keV X-rays. The collaboration has already 
reached an impressive lower limit for the 2i/ECEC decay of T1/2 > 2.6 ■ 10^° y (90%). 
On the other hand, the neutrinoless ECEC decay usually requires an extra photon, 
which must carry away the excess energy given by the Q-value of the transition 
(which is about 2.7 MeV in the ^°^Cd case). This causes an extra suppression and 
the expected half-life is above 10^° years for m^^ = 1 eV by referring to a calculation 
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performed for Sn [10 . 



The ECEC process discussed in this paper brings in a new aspect, which could be 
rather advantageous as far as the experimental conditions are concerned. If Nature 
provides us with a mass degeneracy, as discussed in this paper, the decay proceeds 
through a narrow resonance and the rate increases dramatically giving half-lives 
which are competitive to the neutrinoless /3~/3~ decays. On top of this, there are 
several additional advantages. The de-excitation of the final excited nucleus proceeds 
in most cases through a cascade of easy to detect 7 rays. A two- or even higher- 
fold coincidence setup can cut down any background rate right from the beginning, 
thereby requiring significantly less active or passive shielding. For instance in the 
i36Ba case, the 2315 keV level would go through the cascade 2315(0+) 818(2+) 
g.s.(0"'") (27 rays) and in the ^^^Gd case, the 1952 keV level would de-excite as 
1952(0-) 1242(1-) 89(2+) g.s.(0+) (37 rays). A mere detection of these 7 
rays would already signal the ECEC decay beyond any doubt, as there are no other 
background processes feeding those particular nuclear levels. Further, the lepton 
number conserving ECEC decay with two neutrinos, 

{A, Z) + e,- + e,- -> (A, Z - 2)** + + v,, (4.34) 

is strongly suppressed due to the almost vanishing phase space. Note that the 1v 
phase space has a dependence of ^ with Q being the remaining excess energy. For 



a Q-value in the low MeV region the half-life of the 2z/ECEC is about 10^^ years 42 
For the considered ECEC transitions in fl4.34p the Q-value is assumed to be below a 
keV, which already gives a suppression by at least another 15 orders of magnitude. 
Thus, the signature for a Oz/ECEC resonant transition would not in any way be 
contaminated with 2z/ decays. In fact, the detection of a 7 ray de-excitation cascade 
provides an unambiguous signature for the neutrinoless decay. 



5. Conclusion 

The phenomenon of mixing and oscillations of atoms due to weak interaction 
with the violation of the total lepton number and parity was investigated. We can 
expect that the oscillations can be observed when the energy difference between two 
atomic states is small. The resonant situation can occur, if the daughter atom is in an 
excited atomic or/and nuclear state. The de-excitations of the final state would then 
proceed through the emission of X-rays (atomic structure) and/or 7-rays (nuclear 
structure), both of which could be used as an experimental signature for the process. 
The theoretical framework, which has been employed to describe these processes, is 
similar to that of other oscillation processes, like the oscillations of neutrinos, neutral 
kaons and i?-mesons, or neutrons/antineutrons. 

The concept was applied to the Oz/ECEC process where - because of the final 
state excitation - mother and daughter atoms could be degenerate in mass. An 
experimental detection of such a process would immediately prove the Majorana 
nature of neutrinos. The process can, therefore, be regarded as an alternative to 
the Oz//3~/3~ decay. A list of likely resonant transitions was provided. It was argued 
that accurate mass difference measurements are of paramount importance in order to 
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narrow down the possibilities. In fact, the precision of masses, presently at a several 
keV level, must be improved by at least one or two orders of magnitude. This is 
certainly possible given the progress made over the last years in the use of ion traps. 
Some examples were quoted, where mass measurements have reached accuracies ~ 10 
eV. 

In an ECEC process the daughter atom is left with two electron vacancies, but 
remains otherwise neutral. This extra atomic excitation energy was calculated using 
a relativistic description of the electron wave function, where the Coulomb interaction 
of the two holes was taken into account. In the relativistic description, parity- 
violating transitions from mother to daughter become possible, such as e.g., capture 
of an si/2 and a pi/2 electron leading to the 0"*" — )■ 0+ transition. The electromagnetic 
and Auger decay widths of the atomic states were included in the calculations. In 
the case of mass degeneracy, the decay width determines the rate of transition, and 
the maximum rate corresponds to the unitarity limit. 

The selection rules for electron capture from high orbits were obtained. We have 
shown that nuclear transitions with a change in the nuclear spin A J > 2 are strongly 
suppressed. The capture of K electrons is most likely process due to the large overlap 
of electron wave function with the nucleus. At small distances the wave functions of 
si/2 and pi/2 electrons have similar asymptotes, so that the capture of L2 electrons 
appears to be a feasible process comparable with the K and Li captures. High shells 
are suppressed by a small value of the wave functions in the nucleus, however, these 
states have smaller decay widths, which compensates partly for the smallness of the 
unitary limit of the decay rate of the initial atoms. 

Explicit formulas for the relativistic matrix elements of the Oi/ECEC were given. 

We analyzed the database of the nuclei and their excited states and made a list 
of the most promising candidates for a resonant Oi/ECEC transition. Assuming that 
|m^^| = 1 eV, we arrived at the lower bound of the normalized half-lives ~ 10^^ 
years, which is about a factor of 100 or 1000 below predictions for a time of Oz//3~/3~ 
decay. Because of the uncertainty in the masses, the range of allowed half-lives is 
broad, and reaches several orders of magnitude. Precise measurements of the ground 
state masses, as well as additional spectroscopic information on the excited states 
(energy, spin and parity), are highly warranted to improve predictions of half-lives. 

Our statistical estimates based on the data in Tables |21 El and S] show that 
with the improvement of mass measurements and |wi/3/3| = 1 eV there are about 10 
chances out of 100 to find transition with a normalized half-life below 10^^ years and 
50 chances out of 100 to find transition with a normalized half-life below 10^^ years. 
One can expect two transitions with a normalized half-life below 10^^ years. 

Inverse transitions leading to the neutrinoless production of two bound electrons 
also were considered in cases of a mass degeneracy. The transitions are found to 
be suppressed, because the electrons would have to occupy high orbits, where the 
overlap with the nucleus is small. 

The detection technique for identifying a Oi/ECEC-decay is rather different from 
the Oz//3~/3~-decay, as there is no inherent background from the 2z/-decay, which exist- 
ing experiments have to cope with. Further, by exploiting the coincidence technique, 
in particular if the de-excitation of the nucleus proceeds through a 7-ray cascade. 
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a significantly improved signal to background ratio could be obtained, which would 
alleviate some of the demands on a low-background facility. 

After submission of the manuscript, new experimental results on precision mea- 
surement of Q values and new constraints on the Oi^ECEC half-lives have been re- 



ported [43|, M m |46|, |47 
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Appendix A. Electron wave function inside nucleus 

The relativistic wave function of electron in the Coulomb potential has the form 

where a = (njl), a' = {njl'), /' = 2j — I, and Qjimin) are spherical spinors. 

In the relativistic Coulomb problem, one can distinguish three different scales: 
The first scale is related to the Bohr radius = 1/ (amZ). This scale determines the 
normalization of the electron wave function. The second scale 1/m is the Compton 
wavelength of an electron. Starting from these distances down to zero the upper 
component of Dirac wave function differs markedly from the non-relativistic Coulomb 
wave function. The size of the nucleus is about 500 times smaller than the Compton 
wavelength of an electron, so the effects associated with the finite size of the nucleus 
must be calculated on the basis of the relativistic Dirac equation. One may select 
an even smaller (third) scale. The radius of the nucleus is a few fermi, or turning to 
the proton, its radius is about 1 fm 1/200 MeV~^, which is less than a/m ~ 1/70 
MeV~^. If Z increases, the radius of the nuclei R = 1.2^4^/^ fm increases as well, but 
aZ/m is growing faster, so that R < aZ/m holds for all nuclei. The combination 
of aZ/m defines a third scale below which the Coulomb potential is higher than 
electron mass. 

Asymptotics of the wave functions of electrons work for r < 1/(2A). Note that 
A ~ l/(nas) in the non-relativistic case, and A ~ m in the relativistic case. The 
nuclear radius is below 1 / (2A) and even less than aZ/ (2A). The low r approximation 
should therefore work very well. 

In order to get the electron wave function inside the nucleus we use the asymptotic 
expansion of the Dirac wave function at 2Ar ^ 1: 
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Table A. 5: The upper and lower radial functions of the Dirac bi-spinors, averaged over the volume 
of the nucleus (in keV^^^). 
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2.07x10'' 


2.38x10° 


2.74x10° 


3.48x10° 


4.62x10° 


5.66x10° 


6.31x10° 


<g> 


-5.02x10-2 


-1.75x10-1 


-2.10x10-1 


-2.52x10-1 


-3.46x10-1 


-5.03x10-1 ■ 


-6.49x10-1 


-7.47x10-1 



The normalization and other conventions are those of Ref. [3l|. In particular, 7 = 
y^(J^rT72p"^^^"(aZ)^, A = a/?7?~^^^. The radial functions are real, since aZm — 
kX > 0. Inside the nucleus, the screening of the charge by the electrons is weak, so 
Z is the charge of the unscreened nucleus. 

For uniform distribution of nuclear density, the average values of the upper and 
lower components of the electron wave functions inside the nucleus can easily be 
found: 

3 

< fnjlir) >= ^^-^fnjliR), 

3 

< 9njiir) >= —^gnjiiR). 
At distances r < aZ/m, the ratio (see, e.g., [3l| ) 

MR) _ -2^, j = / + i/2, 

g^j.iR) -f + K^\ f , J = /-l/2, ^ ■ ^ 

is large in the j = 1 + 1/2 states and is small in the j = I — 1/2 states, provided 
aZ <C 1. The electron wave functions with higher values of j are suppressed inside 
the nucleus by additional powers of 2RX <C 1. 

In the non-relativistic theory, wave functions of particles at small distances behave 
like ~ rK Particles, however, become relativistic on the scale of the order of the 
Compton wavelength. As a result, the short distance behavior of the wave function 
changes. Instead of ~ we have ~ r'^~^ (~ r-'"^/^ for aZ <^ 1). The suppression of 
the wave function at short distances is qualitatively different in non-relativistic and 
relativistic theories. 

Table IA.5I shows the upper and lower radial components of the Dirac bi-spinors, 
averaged over the volume of the nucleus. The values are maximum for atoms with 
large Z and for K electrons. One sees a significant suppression of the wave functions 
of 2p3/2 electrons in comparison with lsi/2, 2pi/2, and 2si/2 levels. 

A more accurate calculation of the electron wave functions can be done on the 
basis of relativistic Dirac-Hartree-Fock approximation. 
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Appendix B. Matrix elements for capture of si/2 and pi/2 bound elec- 
trons 



Because of the conservation of angular momentum, transitions from the state O"*" 
to the excited states of nuclei with spin J are possible only when the captured atomic 
electrons have total angular momentum J. The upper and lower components of Dirac 
bi-spinors behave at short distances as ~ r'''~\ where 7 = a/ (j + 1/2)^ — {aZy, j 
is the total angular momentum of electron. Wave functions of electrons with high 
angular momentum are strongly suppressed inside the nucleus. Increase of j by 
one unit leads to suppression of the electron wave function inside the nucleus by an 
amount ~ R/as ^ 1 where R is the size of the nucleus and is the Bohr radius. 
Thus, we restrict ourselves to estimates of matrix elements of neutrinoless double 
capture of si/2 and pi/2 electrons, whose wave functions are given by Eq. f[0]) with 

n,,^(n) = ^ I ^ = J/^' \ = 0' (B.l) 

' ^ ' ^ \ -^o-nxm, J = 1/2, / = 1, ^ ' 

where Xm are the Pauli spinors. We note that the si/2 wave function multiplied by 
75 takes the form of the pi/2 wave function if gg^^^ is replaced by /p^^^ /si/2 
replaced by -g^^i^. 

By taking the advantage of the closure approximation the operator entering the 
nuclear matrix element is of the form of the product of leptonic current and two 
hadronic currents. We have 

E J^(X1)J.(X2) (vl/,^/(xi)C7V(l-75)^^n^,(x2) 

- vl/^^/(xOC7^'^(l - 75)^«m.(x2)) = (1 - (-1)^"+^-^-^) 

X E ^i-. ^.n./(xi)C 7'^7^ (1 - 75)^/.n.,(x2) (xi) (x^) . (B.2) 



m„ra 
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Using the non-relativistic impulse approximation and exploiting the expansion 7/^7^ 
g^lu + i(^ufi we find 



7^J^(xi)7^J''(x2) = -^i^^r„r^5(xi -x„)(5(x2 -x^) x 



n m 



(B.3) 



In what follows, contributions of nuclear matrix elements associated with three terms 
in Eq. (IB. 31) are supplied with subscripts 5* (scalar), tT (time component of the 
tensor), and sT (the spatial component of the tensor), respectively. The nuclear 
matrix elements of Oz/ECEC are presented below. 
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Appendix B.l. 0^^ Oj, 1^ nuclear transitions 

The matrix element of the 0^*" — )■ Oj nuclear transition, entering Eq. f l3.20p . splits 
into three parts 

-M«/3(0+) = MfPiQ^) + Mij)(0+) + Mif ^(0+). (B.4) 

where 

Mlf ^0+) = < 0+ II $^r-r„/.(r„^) [FiJ)(r„,r„) + Gi+)(r„, r„)(r„ ■ r„)] (B.5) 

n m 

X [-4 + (o-n ■ <tJ,] II 0+ >, 

^aP(O^) = ^— < 0+ II V7:-7:„/i(r„„)Gi"t)(r„,r^) - o-^) ■ (r„ x r„) || 0+ >, 

n m 

Mif\0^) = <0+ II ^r-r-/i(r„ J Gi+Hr„,rJ Kx o-J-(r„ xrj II 0+ > . 

n m 

The neutrino exchange potential h{rnm) is defined in Eq. fl3.24p . 

The transition matrix element into the l"*" state has three components 

= M%'\l^) + M^:^\l^) + M^Pil^l (B.6) 

where 

= <l+||5^r-r„Mr„„)Gy(r„,r„)(f„xr„) (B.7) 

n m 

X [-4 + ('^n ■ II > 

Mfp\l^) = z^<l+||5^r-r„M^n™)[Fi-)(r„,r^)(o-„-tT„)-Gi-;(r„,r„) 

X [l^n ■ rm,((T„ — Cr^) + (cr„ — Cr^) ■ r„ f„ — (cr„ — (Trn) ' Tm]] II 0^ > 



n m 



The functions Fj^\rn,rm) and G^^J {rn,rm) depend on the quantum numbers of 
the captured electrons 

I fa{rn)fp{r,n) ± fp{j-n)f^{r-^), (riaSi/2,n/3Si/2), 



\/2 I ±(//3('^n)fl'a('^m) ±fl'a('^n)//3(rm))/2, (n„Si/2, n/3j9l/2) , 

gc{rn)gi3{r„,) ± gi3irn)gairm), 

ga{rn)gp{rm) ± g p{r n) g a{r m) , (?^qSi/2, ri/3Si/2), 
= J i-ga{rn)fis{rm) T f[s{rn)ga{rm))/'2 
V2 I ±(-fi'/3('^n)/a(?^m) =F /a(?^n)fl'/3(?"m))/2, (ri„Si/2 , n^Pl/2) , 

fa{rn)fp{r^) ± fp{rn)fa{rm), (?^aPl/2, ?^/3Pl/2)- 
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If the electrons have the same quantum numbers a = /3, these functions should be 
divided further by 

It is worthwhile to note that for a = ^ (e.g., capture of two K electrons) 
) and Gl^J{rn,rm) vanish. 

Averaging of these functions over the nucleus with the unit weight also gives a zero 
result. However, these functions are multiplied by the nuclear matrix elements that 
are antisymmetric under permutation of the arguments. The outcome, therefore, is 
different from zero, although the naive factorization produces a vanishing result. 

A reasonable estimate can be obtained on the basis of Cauchy's inequality. Let 
aij and bij be antisymmetric tensors. The upper limit of the sum aijbij can be 
evaluated from 



1/2 



1/2 



As an estimate of the matrix element we take the upper limit that splits into a 
product of the lepton and nuclear parts. Combinations of the electron wave functions 
entering the matrix elements of Of — )■ 0^ and 1^ transitions are given in Table 1. 

Appendix B.2. 0^ — OJ, l~j nuclear transitions 

The matrix element decomposes into a sum of two parts that have different tensor 
structure: 



>^a/3(0-) = Mlf (0-) + Mlf)(0- 



where 



Qa 



X 



X 



1 ^n)) 



Cr„, X (Tr, 



(B.8) 
(B.9) 

0+>. 



The transition matrix element of the vector type splits into a sum of three parts: 
M^p{l-) = Mir(l-) + Mf^\l-) + (B.IO) 



where 



M, 



(B.ll) 



X 



— < 17 II ^'^n'^mK^nm) 

n m 



9l 



9a 



X 



Hi/sirn, r^)f„ X - <r„) - H)^p{rm, r„)r^ x (<t^ - <r„) || 0+ >, 



r{-), 



(o-n ■ O-m)] II 0+ >, 
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Mif\l-) = < 17 II }_^r-r„Mr„ J 

n m 



X 



0+> 



The functions ifi^Vrn, r™) have the form 



'a/3 

87rg^^^ (r„,r^) ^ I (^,(r„)^^(r^) ^ /;3(rn)/a(r™))/2 

\/2 I Mgii{rn)ga{rrn)T fc,{rn)fp{rra))l2, (n„Si/2,n^Pi/2), 

[ /a(r„)5'/3(r„) ± //3(r„)5(c,(r^), (^^,^1/2,^/3^1/2). 

If a = /3, then the functions should be divided by \/2. For a = /3 (e.g., capture of two 
electrons from the K-shell), the averaged matrix element of i/|^^^(r„, r^) vanishes. 

It is important to note that due to parity violation, the transitions into the OJ 
states with the capture of two Si/2 electrons are allowed, as well as the transitions 
into the 0^ states with the capture of one Si/2 and one pi/2 electrons. 

The functions if^^''(r„, r^) and if^^''(r„, r^) are explicitly symmetric and anti- 
symmetric with respect to the indices n and m. They are multiplied by the nuclear 
part of the operator, which has the same symmetry properties. Factorization of 
the transition amplitude 0^ — )■ 0~ is possible, since the combination of electron wave 
functions is averaged over the nuclear volume with a unit weight. The nuclear matrix 
element equals f l3.23p . 

In the transition 0^ — I7 we encounter the same problem as in the transition 
0^ — )■ Ij . The nuclear part of the operator is antisymmetric, so the straightforward 
factorization is impossible. We again restrict ourselves to estimate the upper limit 
on the matrix element squared. Before the estimate we must, however, take into ac- 
count that only the antisymmetric part of H^^\rn, r^) contributes to the amphtude, 
therefore, the value of 

(i^y(r„,rJ-/7y(r^,r„))V4 

must be averaged. 

The leptonic parts of the matrix elements discussed above are shown in Table [TJ 
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